Microscopic theory of the two-proton radioactivity 
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We formulate theory of the two-proton radioactivity based on the real-energy continuum shell model. 
This microscopic approach is applied to describe the two-proton decay of the I2 state in ^*Ne. 



o 
o 

00 



> 

in 
o 

o 
\o 
o 



X 



PACS numbers: 23.50.-(-z, 21.60.-n, 24.10.-i 

One of the main frontiers of the nuclear many-body 
problem is the structure of weakly bound and unbound 
nuclei with extreme neutron-to-proton ratios. Weakly 
bound states or resonances cannot be reliably calculated 
in the closed quantum system formalism. For bound 
states, there appears a virtual scattering into the con- 
tinuum phase space which affect the effective nucleon- 
nucleon interaction. For unbound states, the continuum 
structure appears explicitly in the properties of those 
states. Mathematical formulation within the Hilbert 
space of nuclear states embedded in the continuum of 
decay channels goes back to Feshbach A unified de- 
scription of nuclear structure and nuclear reaction as- 
pects is more complicated and became possible in realis- 
tic situations only recently in the framework of the Shell 
Model Embedded in the Continuum (SMEC) j^,^], which 
is based on the completeness of a one-particle basis con- 
sisting of bound orbits and a real-energy continuum. The 
single-particle (s.p.) resonances have to be regularized: 
they are included in a discrete part of the spectrum after 
removing the scattering tails which are fully incorporated 
in the embedding continuum. The configuration mixing 
in the valence space (internal mixing) and through the 
coupling to the scattering continuum (external mixing) is 
calculated microscopically and the asymptotic states are 
obtained in the S'-matrix formalism ^ y| . The SMEC 
contains coupling to one-particle continuum and could 
not be applied applied for the description of Borromean 
systems or the two-nucleon decays. In this Letter, we for- 
mulate the SMEC to include couplings to the two-nucleon 
decay channels and present first application of this new 
formalism to the two-proton (2p) radioactivity. 

Two-particle continuum in SMEC. — Hilbert space is 
divided in three subspaces : Q, P and T. In Q subspace, 
A nucleons are distributed over (quasi-)bound single- 
particle (qbsp) orbits. In P, one nucleon is in the non- 
resonant continuum and A-1 nucleons occupy qbsp or- 
bits. In T, two nucleons are in the non-resonant contin- 
uum and (A-2) are in qbsp orbits. The coupling between 
Q, P and T subspaces changes the 'unperturbed' Shell- 
Model (SM) Hamiltonian (Hqq) in Q into the effective 
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where superscript 
tion, &j^\e) = 
the Green's function in P modified by the coupling to T, 
and G^jt\E) = - Htt]~^ is the Green's function 

in T. In the above equations, Hpp , Htt are the unper- 
turbed Hamiltonians in P, T subspaces, respectively, and 
Hqp, HpQ , HpT, Htp are the corresponding coupling 
terms between Q, P, and T subspaces. The second term 
on the r.h.s. of (Q) describes a di-proton emission, and 
the third term describes a modification due to the mix- 
ing of sequential 2p, di-proton and one-proton (Ip) decay 
modes. In the following, we shall discuss two limits of this 
general process : (i) indirect 2p-emission (sequential 2p- 
emission is a special case of this limit) Htq — Hqt — 0, 
and (it) direct 2p-emission Htp = Hpp = 0. In both lim- 
its, the interference of 2p- and Ip-emissions is taken into 
account by the external mixing of SM wave functions. 

Indirect 2p-emission. — In this limit, the effective 
Hamiltonian ^ can be written in a convenient form 
which separates Ip- and 2p-coupling terms: 



^QQ — EIqq + HqpG^p\e)Hpq 
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HqpG^+\e)Hpt G'f'{E) HtpG'^'{E)Hpq 
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, where {E) = [E^+^ - Hpp] ^ is the Green's func- 
tion in P. To calculate the 2p-emission width of a given 

(A) 

SM state , we begin by diagonalization of first two 
terms on the r.h.s. of (|2J) in the SM basis {'^\'^^}(sm) 

and obtain new many-body states {^'i^^}{ip) in the par- 
ent system which include the configuration mixing due 
to the coupling to the Ip-continuum. From a chosen de- 
caying state '^')'^\ we go on to calculate the width due to 
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the coupling to the 2p-continuum (the last term in JS))). 
In the following, we shall assume that the indirect 2p- 
emission is the sequential process, i.e. the first emitted 
proton is a spectator of the second emission. This im- 
plies a following identification : Hpp — > Hqiqi + p/iop, 
Htt Hp> p> + p/ioPj where primed quantities refer to 
(A-l)-nucleon space, i.e. in Q' subspace, (A-1) nucleons 
are in qbsp orbits and in P' subspace one nucleon is in the 
continuum and (A-2) nucleons are in qbsp orbits, ho is a 
one-body potential describing an average effect of (A-1) 
particles on the emitted proton and p denotes a projector 
on the one-particle continuum states. With this identifi- 
cation, Hpx becomes a coupling between newly defined 
Q' and P' subspaces. 

Direct 2p-emission. — The effective Hamiltonian de- 
scribing a di- proton emission becomes: 



^QQ — HqQ + HqpG^p'{E)Hpq + HQTC'r^' {E)Htq{^) 

As before, we begin by calculating many-body states 
{^\^^}(ip) ill parent nucleus. With the new initial 
state j , we calculate the 2p-decay width: 
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HQTG+{E)HTQ\h.f^)) (4) 



where A^'^^ is an intrinsic state wave function of the par- 



ent nucleus corresponding to ^ j 



\w: 
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the source term, and \ujj''"^') is a continuation of |A^"^'') 
in T. In general, the description of a 2p-decay into T sub- 
space requires a formulation for the three-body asymp- 
totic. In this case, lUj is expanded in hyper-spherical 
harmonics three-body Jacobi coordinate system. The de- 
tails of the CC formalism with three-body asymptotics in 
the SMEC will be published elsewhere • In the follow- 
ing, we shall approximate the 2p-decay by the emission 
of (2p)-cluster and the final-state interaction between the 
two protons in terms of the s-wave phase shift. A similar 
scenario has been used in SM-|-i?-matrix model calcula- 
tions of di-proton decay j5|j|6||. The decay channel is spec- 
ified by : c(C/) = {A^^^"'); [^(fc), (0oo(C/), 5, L)^(^.)]^}, 



where A 



{A-2) 



is the intrinsic state wave function of a 



daughter nucleus corresponding to a non-spurious SM 



J(fc) is the angular momentum of a daugh- 



state <I>^ 

ter system, ^oo(C^) is the Os intrinsic state wave function 
of (2p)-cluster with intrinsic energy U and spin S = 0, L 
is the relative angular momentum of a (2p)-cluster and a 
daughter nucleus, and J is the total angular momentum 
of a system [(A — 2) (g) (2p)]. The source term in is 
expanded in the harmonic oscillator (HO) basis : 
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where r is the relative distance between a daughter nu- 
cleus and a (2p)-cluster and w„^l is the HO wave function. 
The source in ^ is independent of the intrinsic energy 
[/ of a (2p)-cluster. The continuation of A^"^-* in T is 
a solution of the in-homogeneous coupled- channel (CC) 
equations : 



{E - HTT)\i^f^^^) ^ \wj) 
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For Htt, we assume that the total system [(A — 2)® (2p)] 
can be considered as a two-body system in the average 
potential Uq. Projecting ((HJ on the channel c(U) gives : 



E-{Eu + U)-T[^+Uo{r)]T 



T..{+) 
j-AU) 



(r) 



^Ic(C/)(^) 



where E^ is the intrinsic energy of a daughter nu- 
cleus. P^/2/i is the internal kinetic energy of the system 
[(A — 2) (g) (2p)] with the reduced mass /x and T is a pro- 
jection operator on the subspace of cluster states in the 
continuum of P^/2/i -|- C/o(r). The internal energy of the 
cluster U is distributed according to the density of states 
function for the two protons p{U) The di-proton de- 
cay width is then : 
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FIG. 1: Ip-emission width for the decay of 1^ state in ^^Ne 
is plotted as a function of the strength Vo of the WB force 
for the psdfp (solid line) and WBT (dashed line) effective 
interactions. 



Two-proton decay of the 1^ state in ^^Ne. — Nuclear 
decays with three fragments in the final state are very 
exotic processes. The 2p radioactivity is an example of 
such a process which can occur for even-Z nuclei beyond 
the proton drip line : if the sequential decay is ener- 
getically forbidden by pairing correlations, a simultane- 
ous 2p-decay becomes the only possible decay branch. 
The diproton decay may also be observed in a situation 
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where a Ip-decay is allowed, as suggested recently in the 
decay of 1^ state at 6.15 MeV in ^^Ne |3|. Here, the de- 
cay is dominated by a Ip-decay to weakly bound states 
5/2]^ and l/2f in ^^F. No intermediate resonances are 
available in ^^F which would enable a standard indirect 
2p-emission to the ground state (g.s.) of ^^O. The exper- 
imental width for the 2p-emission is 21 ± 3 eV assuming 
a diproton model, and 57 ± 6 eV assuming a sequen- 
tial decay model. In describing the decay of 1^ state in 
"'^^Ne, one has to take into account an influence of the 
Ip-decay process on the 2p-decay through the {5/2+}- 
and {l/2+}-scattering states correlated by weakly bound 
states 5/2+ and 1/2+ of ^^F. 

In SMEC, the radial s.p. wave functions in Q and the 
scattering wave functions in P are generated by a self- 
consistent procedure starting with the average potential 
of Woods- Saxon (WS) type with spin-orbit and Coulomb 
parts included, and taking into account the residual cou- 
pling between Q and P (the procedure of determining 
the self-consistent potentials in the CC equations and the 
regularization of s.p. resonances is described in 0, 
This procedure yields new orthonormalized wave func- 
tions in Q, P and T and new self-consistent potentials 
for each many-body state in Q. In this paper, for the 
effective interaction in Hqq and Hqiqi, we take either 
WBT Hamiltonian Q or USD Hamiltonian for the (sd)- 
shell 0, the KB' interaction for the (p/)-shell ^.nd 
the G matrix for the cross-shell interaction. The 
latter interaction is called the psdfp Hamiltonian. Both 
psdfp- and WBT- Hamiltonians yield the overall correct 
energies of '0?ia;' and 'l?ia;' states in this mass region. 
The residual couplings between Q and the embedding 
continuum is given by the Wigner-Bartlett (WB) force : 
Vi2 = -Vo[a + (1 - a)P[2\S{ri - ra) with a = 0.73 and 
the strength Vq which is adjusted to the Ip-decay width 
of I2 state in ^^Ne. 

The first two terms on the r.h.s. of (0) give an effec- 
tive Hamiltonian describing the Ip-decay and the mix- 
ing of SM states due to the Q - P coupling. The ref- 
erence potential used in the description of Ip-cmission 
of I2 state in ^^Ne ('^^Ne' parametrization) has the ra- 
dius Rq = 3.28 fm, the surface diffuseness a ~ 0.58 
fm, the strength of spin-orbit potential Vgo = 3.68 MeV. 
The depth % = -57.62 MeV of the central part is ad- 
justed to yield the s.p. state 0^5/2 at the experimental 
Ip-separation energy, lpi/2, IP3/2, 0/5/2, 0/7/2 s.p. res- 
onances in this potential are regularized, as described 
in 0, 01 J taking Vcut = 5 fm (6 fm) for 0/ (Ip) reso- 
nances, which corresponds approximately to the top of 
the barrier. The diffuseness of the resonance cutoff is 1 
fm. This reference WS potential 0,0] defines energies of 
s.p. states and yields radial s.p. wave functions for those 
channels which are not affected by the continuum cou- 
pling. For all other channels, the diagonal potentials in 
CC equations which include the continuum-coupling cor- 
rection are constrained to reproduce reference potential 



s.p. energies ^ig. 1 exhibits the Ip decay width 

of I2 state as a function of the strength Vq of the WB 
coupling to the embedding continuum. The Ip width 
saturates for large Vq which is a well-known feature of 
open quantum systems. The calculated width is smaller 
than the experimental value 50 ±5 keV. In the following, 
we shall take Vq = 900 MeV-fm^ for the strength of the 
coupling between Q, P and Q, T subspaces. 
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FIG. 2: The real part of a di-proton source function for the 
decay of 1^ state in ^®Ne. Solid (dashed) line corresponds to 
the psdfp (WBT) effective interaction. The short-dashed line 
denotes the source function obtained by neglecting couplings 
to opened Ip-emission channels. For more details, see the 
description in the text. 

The sequential 2p-emission to the g.s. of ^^O is possi- 
ble through the correlated continuum states of ^^F. The 
theoretical scheme is described above. In a description 
of the first proton emission into 5/2+ and 1/2+ contin- 
uum of ^^F, we use the '^'^Nc'-refcrcncc potential. Refer- 
ence WS potential used in the description of the emission 
of the second proton ('^^F' parametrization) has the ra- 
dius Rq — 3.21 fm and the surface diffuseness a — 0.58 
fm. The strength of the spin-orbit potential (V^o = 3.68 
MeV) and the depth of the central part {Vq = -52.46) 
MeV are adjusted to reproduce the experimental Ip- 
separation energies in ^''F for s.p. states 0c?5/2, lsi/2. Po- 
tential Hq describing an average effect of (A-1) particles 
on the first emitted proton is given by the self-consistent 
potential obtained from '^^Ne' reference potential. 

The width for sequential 2p-emission r|2p)'' depends 
strongly both on the effective SM interaction and on the 
coupling to Ip-decay channels. If the coupling to the 
Ip-scattering continuum is neglected, i.e. the external 
mixing of SM states is neglected, then : Plgp)"* = 13-1 
and 38 eV for psdfp and WBT interactions, respectively. 
The branching ratio B1/2+ for the 2p-decay through the 
{l/2+}-continuum in ^'F is 87% and 96% for these two 
interactions, rf^*^?-* and B1/2+ depend on the total in- 
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tensity of {lsQd){lpOf) component in 1^ state which is 
5.8% and 11.3% for psdfp and WBT interactions, re- 
spectively. Including couphngs to the Ip-decay channels 
one finds : r|2p'^ = 88.8 and 13.6 cV for psdfp and 
WBT interactions, respectively, i.e. a bigger 2p width 
is obtained using a psdfp interaction for which the in- 
tensity of (Is0d)(lp0/) component in the 1^ SM state is 
smaller than for the WBT interaction but the external 
mixing is much stronger than for the WBT interaction 
and, moreover, interferes constructively. The intensity 
of the {Is0d){lp0f) component in ^2 (ip) SMEC state is 
now 18.9% for psdfp interaction and 8.7% for WBT inter- 
action. Similarly, the external coupling is reversing the 
tendency for branching ratio -B1/2+ which becomes 93% 
and 80% for psdfp and WBT interactions, respectively. 
Clearly, the coupling to opened Ip-decay channels is es- 
sential for understanding the sequential 2p-decay of I2 
state in ^^Ne. 
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teraction neglecting coupling to Ip-decay channels. All 
three source functions are qualitatively similar and the 
Q - P coupling strongly modifies the magnitude of the 
source. The imaginary part of the source (see Fig. |2Jl is 
due to the coupling to Ip-emission channels and strongly 
differs for psdfp and WBT interactions. The width for a 
direct 2p-emission r^p is 1.89 and 1.01 eV for psdfp and 
WBT interactions, respectively. Neglecting the coupling 
to Ip-decay channels, one finds 0.8 and 1.17 eV. Similarly 
as for the sequential 2p-decay, Tj*'^'' depends strongly on 
the coupling to the opened Ip-decay channels. 

In summary, we have extended the real-energy con- 
tinuum shell model (SMEC) to include couphngs to the 
two-particle continuum and applied this new formalism 
for Ip- and 2p-decays of 1^ state in ^*Ne. The ex- 
perimental data are compatible with the sequential 2p- 
decay through the correlated scattering continuum of ^ ''F 
with, possibly, a weak di-proton branch. As compared 
to the SM-t-R-matrix calculations , we find systemati- 
cally larger (resp. smaller) two-proton decay width 

(resp. Fjp"^^). This difference is mainly due to the ab- 
sence of external mixing in SM-|-R-matrix model |^ 0| 
and a more realistic description of the emission process 
in SMEC. The 2p-decay width is strongly influenced by 
the interference between external and internal mixings of 
SM wave functions. Strong couplings between Ip- and 
2p- emission sectors invalidate a simple picture of direct 
2p-decay as a new independent decay mode, at least if 
the Ip-decay channels are opened. 

We wish to thank F.C. Barker and F. Nowacki for use- 
ful discussions. 



FIG. 3: The imaginary part of a di-proton source function 
for the decay of 1^ state in ^®Ne. For more details, see the 
caption of Fig. Inland the description in the text. 

The theoretical scheme used for describing a di-proton 
decay is given above. Radial s.p. wave functions which 
enter in the calculation of the source function wj^{r) are 
generated using '^^Ne'-reference potential. In the studied 
case, the (2p)-cluster is emitted with a relative angular 
momentum L = 1. To calculate ujj we suppose that 
the interaction of the (2p)-cluster with a daughter nu- 
cleus is described by an average potential Uq which 
is a sum of central WS and Coulomb potentials. Pa- 
rameters of Uq are deduced from the deuteron scattering 
data 113. The depth of the WS part of Uq {% = -57.97 
MeV) is adjusted to obtain a p-wave resonance for a par- 
ticle of mass 2TOp and charge Z = 2 at the energy avail- 
able for the 2p-decay of 1^ state in ^^Ne. Fig. |2] shows 
the real part of a di-proton source which for psdfp (the 
solid line) and WBT (the dashed line) interactions. The 
short-dashed fine in Fig. |21 shows results for a psdfp in- 
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